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INTRODUCTION 

Shaped  charge  jets  are  being  used  as  initiating  stimuli  in  tests  of  energetic 
materials.  The  response  of  these  materials  can  vary  considerably,  depending  on 
which  shaped  charge  is  used  in  the  test.  One  jet  parameter  which  seems  to  be  of 
importance  is  the  energy  that  the  jet  deposits  in  the  target.  This  energy  is  spread  out 
over  the  penetration  path  in  a  complicated  way  and  a  means  to  calculate  this  energy 
profile  was  needed.  Thus  the  hydrodynamic  penetration  theory  was  employed  to 
calculate  these  energy  deposition  profiles  to  assess  effects  of  standoff  and  various 
combinations  of  conditioning  plates  on  the  stimulus  level  seen  by  test  samples. 

METHOD 


Calculations  of  energy  deposition  were  carried  out  using  the  hydrodynamic 
penetration  model.  The  energy  deposited  by  the  jet  in  penetrating  an  infinitesimal 
distance  dp  is  just  the  kinetic  energy  of  the  jet  which  is  used  up  in  obtaining  this  dp. 
Thus,  the  energy  deposition  can  be  shown  to  be 


iiE= 

dp 


(1) 


where  Uj  and  Aj  are  instantaneous  values  of  jet  velocity  and  cross-sectional  area  at  the 
stagnation  point  of  the  penetration  cavity  while  h  and  r,  are  the  respective  jet  and  target 
densities. 


The  cross-sectionai  area  of  the  jet  constantly  changes  due  to  jet  stretching.  At  a 
fixed  Lagrangian  point  on  the  jet  (which  we  choose  to  measure  using  the  accumulated 
jet  mass),  the  area  change  is  given  by 


rjAj 


2dUj 

dm 


where  m  is  the  accumulated  mass  in  the  jet,  measured  back  from  the  tip.  If  we  invoke 
the  usual  assumption  that  the  velocity  of  a  jet  particle  does  not  change  during  flight, 
then  dU/dm  i.s  not  a  function  of  t  and  the  jet  area  is  given  by 


J  dU. 

(2) 

with  integration  constants  Aq  and  to .  While  penetrating  a  target,  the  rate  at  which  mass 
is  used  up  is  given  by 


llia=r.A.— 
dt  J  Jl+p  J 


(3) 


1 


with 


After  subontution  of  equation  S.'the  expression  for  mass  loss  Is 


dm,, 

dt 


P 

1  +  p 


(4) 


where  Uj ,  Ao  and  dUj/dm  are  only  functions  of  m.  A. simple  integral  of  penetration 
velocity  gives  penetration  depth, 


dt 


(5) 


In  this  formulation,  the  shaped  charge  jet  is  completely  described  by  stipulating 
velocity  and  cross-sectional  area  at  a  given  time  as  a  function  of  arcumulated  mass. 
There  are  numerous  ways  to  handle  the  required  initial  cpndjjions.  One  possibility  is 
to  calculate  the^cross  section  when  the  element  of  jet'fwms  and  Record  the  time  when  if 
forms  as  a  function  of  accumulated  triass  tip  to  this  point.  Another  approacH  is  to  use.a 
jet  snapshot  to  measure  area,  in  wfij^  c£^  iojis  fhe  same  for  eacli  jet  eiemeht.  This 
latter  method  was  used  here  and  the  time  of  the  snk^fiot  Ts  relative  to  the  instant 
when  the  firstjet  particle  forms.  ’  '  ’ 


Note  that  the  equation  2  can  be  used  to  calculate  (dU)/dm)  if..tw6  x-ray 
snapshots  of  a  jet  are  taken  at  different  times.  The  time  inteival  is  chosen  large 
enough  that  the  jet  area  change  is  easy  to  measure.  With  area  from  the  earliest 
snapshot  and  (dUj/dm)  from  equation  2,  a  cbipplete  initial  condition  data  set  is 
established. 


VIRTUAL  ORIGIN  ASSUMPTION 


Some  simplifications  can  be  obtained  if  the  virtual  origin  jet  assumption  is 
invoked.  A  virtual  origin  jet  has  a  linear  velocity  distribution  so  that 


3U. 

or  U.=-f(t)x  +  U,p 

At  any  instant  of  time,  to 


(6) 


9m  3ra  9x  .  /  i  \ 
9U.  “  9x  9U.  "'^j^o('fi^j 


2 


However,  accumulated  mass  is  a  Lagrangian  coordfnate  so  that 


3m  dm 
3u]“dUj 


and  thus 

.  dUj  ^to) 

is  a  constant.  Then  equation  4  can  be  integrated  to  give  the  velocity  at  the  bottom  of 
the  penetration  cavity,  viz. 


^=[1 

Now  equation  5  is  integrated  using  this  expression  to  give 


Choose  to  as  the  time.when  the  jet  tip  hits  the  target.  The  time  from  the  virtuai  origin  to 
the  target  is  found  from 


*  _  ^lip“ 


with  (')r  referring  to  the.last  element  of  the  jet,  which,  after  substituting  for  the  velocity 
difference  from  equation  6  leads  to 


giving  the  usual  form  of  the  Allison  -  Vital!  formula  (ref  1,  p.138),  viz. 


P  =  Uupto 


t  w 


-1 
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Calculations  of  Energy  Deposition  Rates  for  Several  Jets 

Ernest  Baker  of  ARDEC  provided  a  copy  of  the  analytical  shaped  charge  cods 
PASCC1  (ref  2),  which  he  developed,  (and  also  tutorials  to  the  author  which  allowed 
successful  use  of  his  code)  to  calculate  jet  formation  conditions  for  several  changes. 
Another  code  was  written  to  integrate  equations  4  and  5  using  the  output  from 
PASCC1.  Some  results  are  discussed  belovr. 

First,  a  viper  jet  was  calculated  using  PASCC1.  Insufficient  detail  of  the  cone 
apex  was  included  for  PASCC1  to  predict  any  tip  particle  phenomena  for  this  jet.  The 
calculated  effect  of  standoff  on  energy  dissipation  for  a  Viper  shaped  charge  against  a 
homogeneous  target  with  density  of  1.1  gm/cc  located  at  36.83  cm  (14.5  inches)  from 
the  base  of  the  Viper  is  shovm  in  figure  1.  A  curious  feature  is  that  energy  deposition 
peaks  at  some  distance  into  the  target  because  jet  area  increases  at  the  stagnation 
point  more  rapidly  than  the  square  of  velocity  falls  off,  according  to  the  calculation. 
Standoff  has  almost  no  effect  after  about  25  cm  of  target  has  been  penetrated. 

The  calculated  effect  of  various  thickness  of  steel  conditioning  plate  is  shown  in 
figure  2.  These  plates  are  assumed  to  be  located  12.7  cm  (five  inches)  from  the 
charge,  while  the  1 .1  gm/cc  target  is  still  located  at  36.83  cm.  Thin  plates  wipe  off  the 
early  portion  of  the  jet  causing  the  large  area  portion  to  hit  target  material  at  a  lesser 
distance  than  i!  the  early  material  were  allowed  to  impact  the  target.  The  jet  cross 
section  has  less  time  to  shrink,  which  causes  the  maximum  value  of  energy  dissipation 
to  be  higher  with  the  thinest  plate  than  without  it.  It  becomes  clear  that  energy 
deposition  is  sensitive  to  details  of  the  jet  structure. 

A  comparison  of  three  typical  charges,  all  at  36.83  cm  standoff  from  the  1.1 
grn/cc  target  is  shown  in  figure  3.  Energy  deposition  by  the  BRL-81  mm  runs  about 
50%  above  the  Viper.  The  Tow2a  produces  a  jet  with  wildly  varying  radius.  The  thin 
early  jet  cannot  deliver  much  energy  to  the  target,  but  the  rate  of  energy  delivery 
Increases  rapidly.  There  are  large  particles  in  the  Tow2a  jet,  caused  by  reverse 
velocity  gradients  created  during  jet  formation,  which  travel  without  stretching  and 
deliver  large  energy  input  to  the  target  as  they  impact.  The  figure  shows  that  two  of 
these  particles  should  impact  within  the  0.5  meters  thickness  of  the  target.  Obviously, 
the  Tow2a  delivers  a  completely  different  energy  imput  profile  as  compared  to  the 
other  tv/o  charges  and  the  particles  dominate  energy  deposition. 

Siccit  the  virtual  origin  approximation  is  so  widely  used,  it  is  worthwhile  to 
calcuiate.  tlur  energy  deposited  by  a  virtual  origin  jet  having  the  Viper's  velocity 
gradient  and  mass. '  in  order  to  calculate  energy  deposition  for  virtual  origin  jets,  a 
knowlege  cf  how  their  cross-sectional  area  varies  along  the  jet  is  needed.  As  an 


4 


example,  pick  the  simplest  possible  case,  i.e.,  assume  that  the  jet  has  constant  radius 
along  its  length  at  any  instant  of  time,  then  its  energy  deposition  is  easily  calculated, 
and  is  given  by 


The  value  of  AoS  is  calculated  from 


(7) 


1  U 


m 


v/ith  m  the  accumulated  mass  up  to  the  particular  element  chosen  as  the  rearmost 
element  (note  that  this  product  is  constant  as  the  jet  stretches).  The  standoff,  S,  is 
calculated  from  the  PASCC1  Viper  output  in  the  following  way.  The  PASCC1  output 
gives  jet  velocity  as  a  function  of  axial  distance  as  a  snapshot  picture  at  some 
particuiar  time.  Let  xi  be  the  distance  of  the  tip  from  the  cone  reference  point,  the  point 
where  the  first  jet  particie  is  formed.  The  distance  of  the  tip  from  its  virtual  origin  can  be 
found  from 


where  h  is  the  distance  from  the  tip  to  the  chosen  rear  element  which  has  velocity  Ur 
(taken  nominally  as  2.5  K/S ).  Then,  if  xa  is  the  distance  from  charge  base  to  the  target 
and  X3  is  the  distance  from  the  charge  base  to  the  cone  reference  point,  the  target 
standoff  distance  is 


S  =  Xj  +  Xj  +S  -  Xj 


Formula  7  is  plotted  in  figure  4  aiong  wrtth  the  calculated  results  from  PASCC1 .  The 
constant  area,  virtual  origin  jet  deviates  from  the  PASCC1  jet  for  the  first  15  cm  of 
target  penetration  and  this  deviation  is  considerabie.  Net  deposited  energy  generally 
varies  with  the  square  of  tip  veiocity,  but  is  notably  sensitive  to  the  area  distribution 
along  the  jet. 

Experimental  Viper  Jet 


Mr.  James  Pham  of  ARDEC  provided  a  double  exposed  x-ray  of  a  viper  jet  for 
analysis.  The  x-ray  shows  the  jet  at  28.904  microseconds  and  49.329  microseconds 
reialive  to  the  instant  that  the  first  piece  of  jet  forms.  Note  that  the  firing  pulse  occurs  at 


-12.146  microseconds  on  this  time  scale,  so  the  time  relative  to  the  fitinci  pulse  has  this 
value  added  on.  Both  jet  images  were  digitized  and  the  results  are  shown  in  figures  5 
and  6.  A  best  estimate  for  the  axis  of  symmetry  was  made  by  fitting  curves  between 
the  points  away  from  the  tip.  As  can  be  seen,  the  resulting  axis  nicely  pass  through 
the  center  of  the  tip  also.  Based  on  these  axis,  the  axisymmetric  volume  at  each  time 
was  determined  as  a  function  of  distance  from  the  tip.  The  resulting  jet  radius  and  jet 
mass  (volume  times  density)  are  shown  in  figures  7  and  8.  Consistent  with  the 
hydrodynamic  development  used  here,  it  is  assumed  that  the  point  described  by  "m" 
grams  back  from  the  tip  of  one  X’  ray  is  also  "m“  grams  back  in  the  other.  Thus,  the 
assumption  is  that  no  inversion  of  mass  coordinates  occurs  during  flight.  The 
difference  in  radii  between  x-rays  is  nominally  20%,  so  that  small  errors  in  reading  the 
data  are  not  crucial.  Equation  2  can  now  be  used  to  calculate  dUj/dm  as  a  function  of 
accumulated  mass  (fig.  9)  and  ultimately  jet  velocity  (fig.  10).  The  range  of  the 
accumulated  mass  coordinate  is  limited  by  how  much  of  the  early  jet  appears  in  the  x- 
ray.  Results  of  this  procedure  give  a  smooth  velocity  profile  with  two  small  regions  of 
reverse  gradient.  The  x-ray  tip  velocity  is  faster  than  that  calculated  by  PASCC1  as 
can  be  seen  when  comparing  velocity  distributions  in  figure  1 1 .  The  velocity  gradient 
after  mass  coordinate  =.002  kg  is  similar  to,  but  slightly  steeper  than  that  calculated  by 
PASCC1 .  A  comparison  of  the  jot  radii  from  the  x-ray  with  that  calculated  by  PASCC1 
is  shown  in  figure  12.  The  x-ray  results  appear  to  agree  with  the  calculation  except  in 
tlie  tip  region.  Energy  deposition,  calculated  using  these  experimental  data,  is  shown 
on  figure  13.  Obviously,  the. tip  delivers  a  lot  of  energy.  Energy  deposition  drops  to 
that  from  PASCC1  after  the  tip  is  absorbed.  These  results  on  a  larger  scale  are 
shown  in  figure  14.  It  seems  the  PASCC1  resuits  can  be  used  satisfactorily  beyond 
about  4  cm  penetration.  The  effect  of  a  1.27  cm  steel  conditioning  plate  on  the  x-ray 
derived  jet  data  is  shown  in  figure  15.  The  large  lip  effect  is  essentially  removed  by  this 
conditioning  plate. 

CONCLUSION 

With  the  use  of  various  shaped  charges  as  initiating  stimuli  for  probing  assorted 
energetic  materials,  the  ability  to  quantify  the  strength  of  stimulus  becomes  important 
for  comparing  results  and  in  forming  a  basis  for  understanding  what  parameters 
govern  the  response  of  the  energetic  material.  This  hydrodynamic  analysis  gives  a 
quick  way  to  quantify  the  stimulus  for  any  test  setup  and  also  serves  to  point  up  how 
sensitive  the  stimulus  parameters  can  be  to  details  of  the  jet.  It  is  particularly  useful  to 
study  cases  where  stimulated  reaction  appears  only  along  part  of  the  path  of  the  jet 
through  the  target.  The  point  of  cessation  of  reaction  can  be  correlated  with  local  jet 
properties.  Of  course,  this  analysis  is  only  as  good  as  the  basic  jet  description  and 
developing  an  accurate  jet  description  should  be  pursued. 


Figure  2.  Effect  of  conditioning  plates 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

CO 

(O 

to 

CO 

CM 

(l/\I/P)  ABjeua 


10 


Figure  4.  Viper  comparisons 
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RgureS.  Jet  at  41.05  microseconds 
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Figure  6.  Jet  at  61 .475  microseconds 
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Figurs  7.  Analysis  of  jet  at  41 .05  microseconds 


Figure  8.  Analysis  of  jet  at  61.475  microseconds 
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Figure  10.  Particle  velocity 
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Viper  velocity  distribution 


♦ 


(W)  sniissi 


18 


h  OOO'O 


^oerat  14.5  inches 


20 


Rgurs  14.  Viper  comparisons 


Viper  with  1/2  inch  plate 
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